Searching for modifications to the exponential radioactive decay law with the Cassini 

spacecraft 
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Data from the power output of the radioisotope thermoelectric generators aboard the Cassini 
spacecraft are used to test the conjecture that small deviations observed in terrestrial measurements 
of the exponential radioactive decay law are correlated with the Earth-Sun distance. No significant 
deviations from exponential decay are observed over a range of 0.7 — 1.6A.U. A 90% CI upper limit 
of 0.84 X 10"* is set on a term in the decay rate of ^^^Pu proportional to 1/E? and 0.99 x 10 * for 
a term proportional to 1/R. 

PACS numbers: 23.60.+e, 23.40.-s, 95.55.Pe, 96.60.Vg, 0620. Jr. 



A recent archive preprint reports evidence for a cor- 
relation between nuclear decay rates and the Earth-Sun 
distance[H. This correlation is extracted from an an- 
nual modulation in the observed decay rates oi^^ Si/^^Cl, 
both (3 emitters, and ^^^Ra, an a emitter. Reference 1 
analyzes this as a correlation with 1/R^{t), the variation 
in the Earth-Sun distance due to the eccentricity of the 
Earth's orbit. To set a scale for this correlation the am- 
plitude of the decay rate variation is ~ 0.1% roughly in 
phase with the 3% annual modulation in 1/R^{t) sug- 
gesting a (3 X 10~^)/i?^ term in the decay rate. 

In the conclusion of this paper the authors observe 
that: These conclusions can be tested... [by] measure- 
ments on radioactive samples carried aboard spacecraft 
to other planets [which] would be very useful since the 
sample-Sun distance would vary over a much wider range. 
I report here the results of exactly such a measurement 
based on the power output of the Radioisotope Thermo- 
electric Generators (RTG) aboard the Cassini spacecraft 
which launched in 1997 and reached Saturn in 2004. 

Cassini is powered by three RTGs each of which is a 
very large (7.7Kg, ISOKCu) "^^^Pu radioactive source, 
an a emitter with an 87. 7y half life 3- The heat from 
these sources are converted to electric power with ther- 
moelectric piles. Together these sources produced 878w 
of electrical power from ^ ISKw of radioactive decay 
heat at launch. The power output of these RTGs were 
literally the lifeblood of the Cassini mission. Their power 
output was monitored carefully and often. 

The trajectory of the Cassini spacecraft is available 
on the web[l]. Pve used these data, converted to astro- 
nomical units {A.U.), to compute l/i?(i)^ = i?g/[x(t)^ -I- 
yit)^ + z(<)2]. Over the first 2 years Cassini went from 
R(0) = 1 at launch, made 2 visits to Venus at R=0.7 and 
crossed the orbit of the earth a total of 4 times before 
finally gaining enough speed to reach Saturn at R=9. 

JPL kindly provided^] P(t), the total electrical power 
from the three RTGs aboard, measured daily since launch 
and the expected power output from their RTG mod- 
eling. The distance of Cassini from the Sun and the 
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FIG. 1: Left, suppressed zero, scale: solid red curve. Helio- 
centric distance [R{t){A.U.)], Black points are the 5 times 
when R(t) = l (green line). Right, suppressed zero, scale: blue 
diamonds, RTG electrical power. 



electrical power output are plotted in Figure 1 for the 
first 2 years of the mission. R ranged from 0.7 — 1.6A.U. 
(0.35 < l/i?2 < 2.2). The power dropped from 878w at 
launch to 815w over this period. 

The thermal power output of an RTG is directly pro- 
portional to the decay rate of the radioisotope generating 
the heat: Pth{t) — No\exp{~Xt)Ed, where Ed is the en- 
ergy released per decay. The electrical power output, 
P{t) = Pth{t)eoe{t) is modified by an initial efficiency, eo, 
and a time (or power) dependent thermoelectric conver- 
sion efficiency; e(i),e(0) = 1. 

I cannot safely use the RTG efficiency model here 
since it assumed only exponential behavior for radioac- 
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FIG. 2: RTG power; green open circles, decay corrected 
[e{t)U{t)], black solid circles, efficiency and decay corrected 
[C/(t)], cyan curve, 3rd order polynomial correction function 
described the in text, red solid line, expected eflFect extrapo- 
lated from Reference 1: 1 - {l/R^{t) - l)/30. 
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FIG. 3: Black points; A{t), RTG normalized thermal power 
less a 3rd order polynomial correction. Red curve; 10 times 
the 90% CL limit, lOa/7?^, Blue curve; 10 times the 90% CL 
limit, 10/3/-R. 



tive decay. Any new physics effect might be inadvertently 
subsumed into that model. The Cassini trajectory pro- 
vides a natural calibration for e(i) using the measured 
power at the 5 points (shown in Figure 1) where the 
spacecraft was lA.U. from the sun. Fitting these mea- 
surements to PR=i{t) = p[0)2(-*/^''-'^v)2(-t/T.jf) yields 
Teff = 21.2 ± 1.9?/. This simple model agrees in shape 
with Reference 4 and obeys the requirements of Carnot 
thermal efficiency: as the power decreases, and the tem- 
perature difference across the thermoelectric piles de- 
crease, the efficiency can only decrease. Something as 
complicated as a space-born thermoelectric pile requires 
more than one parameter to accurately describe its be- 
havior. The single exponential reduces a 5% power drop 
in the first 2 years to a ±1% variation from unity. 

Following Reference 1 I've plotted e{t)U{t) = 
P {t)/P {0)2^ /^'^■'^y, the normalized electrical power cor- 
rected for "^^^Pu decay, and U(t), the normalized thermal 
power, as a functions of time since launch in Figure 2 for 
the first 2 years of Cassini's voyage to Saturn. 

To compare to the rough magnitude of the effect 
reported by reference 1 Fve plotted Urefi{t) = 1-1- 
(0.1%/3%)(l/i?(t)2 - 1) on Figure 2 to extrapolate the 
small R variation of reference 1 for comparison with the 
larger R range available in this study. Extrapolating a 
0.1% decay rate change for a 3% change in to a 50% 
change in 1/i?^ {R = 0.7 A.U., t = 0.43y) should cause a 



-1-4% change in the power output of the RTG 5 months 
after launch. In fact U{t) decreases by ~ 1%. Changes 
this large are excluded by the Cassini efficiency corrected 
data both in magnitude and by the absence of any refiec- 
tion of the shape of 1/R^{t) in either normalized power 
curve. 

In order to set quantitative limits I've fit the efficiency 
and half-life corrected U(t) normalized power data (sup- 
pressing the two obvious outlying points in Figure 1) to 
Psit), a 3th order polynomial in time, to phenomeno- 
logically describe the last 1% variation in the U{t). As 
shown in Figure 2, this polynomial smoothly interpolates 
the U{t) measurements. A 3th order polynomial make a 
very poor fit to 1/R'^{t); these two shapes are approxi- 
mately orthogonal. The error assigned to each U(t) mea- 
surement by requiring this fit to have x^/'^ = 1 is 0.0015. 
The individual relative power measurements have a res- 
olution of 0.15% and the polynomial is a ~ 6ct system- 
atic correction beyond the simple exponential efficiency 
model. 

The difference of the data from the polynomial fit; 
A{t) = U{t) — Psit) are plotted in Figure 3. Some struc- 
ture at the 1(7 level and some outlying measurements 
remain. These difference data are fit to a/R'^ and (3/R 
to give limits on the contribution of a term in the ^^^Pu 
decay rate dependent on the Earth-Sun distance. The 
90%CL fimit on from these fits are \a\ < 0.84 x 10"'' and 
/3| < 0.99 X 10^** respectively. The limiting functions, 
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scale up by a factor of 10 for visibility, are also shown in 
Figure 3. a is to be compared with the correlation seen 
in Reference 1 for ^^^Rn decay of ~ +3 x 10~^. 

The Cassini RTG power data exclude any variation of 
the "^^^Pu nuclear decay rate correlated with the distance 
of the source from the Sun to a level 350 x smaller than 
the effect reported by Reference 1. ^^^Pu and ^^^Ra 
are similar a emitters. Another physical or experimental 
cause of the reported annual variations in nuclear decay 
rates appears to be necessary. More generally Ruther- 
ford, Chadwick, and Ellis's 1930 conclusion that The rate 
of transformation of an element has been found to be con- 
stant under all conditions.^ now has solid experimental 
support at least from Venus (R=0.7) to Mars (R=1.5). 

I am indebted to several of my colleagues for calling 
this paper and physics issue to my attention; Chris Quigg 
and Martin Hu of Fermilab and Jurgen Engelfried of the 
Universidad Aut0noma de San Luis Potosf, Mexico. I am 
also thankful for very useful lunchtime conversations on 
this subject with several of my Fermilab colleagues. I an 
indebted to Richard Ewell and Torrence Johnson of JPL 



for making the RTG data available and Stephen Parke of 
Fermilab for critical comments on this manuscript. 
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